SLAC-PUB-7258 



hcp-ph/9609214 



August 1996 



LIFETIME DIFFERENCE OF Bs MESONS AND ITS 
IMPLICATIONS " 



M. BENEKE 

Stanford Linear Accelerator Center, 
Stanford University, Stanford, CA 94309, U.S.A. 

We discuss the calculation of the width difference AF^^ between the Bs mass eigen- 
states to next-to-leading order in the heavy quark expansion, l/m^-corrections are 
estimated to reduce the leading order result by typically 30%. The error of the 
present estimate (Ar/r)^^ = 0.16^Q Qg could be substantially improved by pin- 
ning down the value of {BB\(biSi)g_p(bjSj)g_p\Bs) and an accuracy of 10% in 
{Ar/r)g^ should eventually be reached. We briefly mention strategies to measure 
{Ar/r)fl^, and its implications for constraints on AMb^, CKM parameters and 
the observation of CP violation in untagged Bs samples. 



To appear in the Proceedings of the 
28th International Conference 
on High Energy Physics, 
Warsaw, Poland, 25-31 July 1996 



"Research supported by the Department of Energy under contract DE-AC03-76SF00515. 



1 



1 Introduction 



Mixing phenomena in neutral B meson systems provide us with an impor- 
tant probe of standard model flavordynamics and its interplay with the strong 
interaction. As is well-known, non-zero off-diagonal elements of the mixing 
matrix in the flavor basis {\Bs), \Bs)} are generated in second order in the 
weak interaction through 'box diagrams'. In the Bg system^ the off-diagonal 
elements obey the pattern 

The mass and lifetime difference between eigenstates are given by ('H' for 
'heavy', 'L' for 'light') 

AMb^ = Mh - A/l = 2 I Afi2 1 , (2) 
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ATb^=Tl-Vh- pg-^ 2ri2, (3) 

up to very small corrections (assuming standard model CP violation). Antic- 
ipating the magnitudes of the eigenvalues, we have defined both AMb^ and 
AFb^ to be positive. Note that the lighter state is CP even and decays more 
rapidly than the heavier state. 

The lifetime difference is an interesting quantity in several respects. Con- 
trary to the neutral kaon system, it is calculable by short-distance methods 
and directly probes the spectator quark dynamics which generates lifetime dif- 
ferences among all b hadrons. If the mass difference AMb^ turns out to be 
large, the lifetime difference also tends to be large and may well be the first 
direct observation of mixing for Bs mesons. If ATb^ is sizeable, CP violation 
in the Bs system can be observed without flavor-taggingOJ. _ 

The following sections summarize the calculation of Ref.El and discuss some 
of the implications of a non-zero AF b^ ■ 



2 Heavy quark expansion of AFs^ 

The mass difference is dominated by the top-quark box diagram, which reduces 
to a local AB = 2 vertex on a momentum scale smaller than M^y . The lifetime 
difference, on the other hand, is generated by real intermediate states and is 

''For Bii mesons there is further CKM suppression and their lifetime difference will not 
be considered here. 
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not yet local on this scale. But the b quark mass mf, provides an additional 
short-distance scale that leads to a large energy release (compared to Aqcd) 
into the intermediate states. Thus, at typical hadronic scales the decay is 
again a local process. The lifetime difference can then be treated by the same 
operator product expansion that applies to the average Bg lifetime and other 
b hadronscl. 

Summing over all intermediate states, the off-diagonal element of the 
decay width matrix is given by 

= 2A^^^^' ^""^ jd^^THeff{x)HeffmBs) (4) 



with 



V2 

+ C2{^l){biCi)v-A{cjSJ)v-A)■ (5) 

Cabibbo suppressed and penguin operators in TLeff have not been written 
explicitly. In leading logarithmic approximation, the Wilson coefficients are 
given by 6*2,1 = (C+ ± C_)/2, where 



as{Mw) 
as(p) 



6/23 



-12/23 



(6) 



and /i is of order to^ . 

The heavy quark expansion expresses AFs^ as a series in local AB = 2- 
operators. In the following we keep l/wfj-corrections to the leading term in the 
expansion. Keeping these terms fixes various ambiguities of the leading order 
calculation, such as whether the quark mass mb or meson mass Mb^ should be 
used, Bjad establishes the reliability of the leading order expression obtained 
in Ref.tttl. Compared to the 'exclusive approach' pursued in Ref.B that adds 
the contributions to AF^, from individual intermediate states, the inclusive 
approach is model-independent. The operator product expansion provides a 
systematic approximation in Aqq d /™bj but it relies on the assumption of 'local 
duality'. The accuracy to which one should expect duality to hold is difficult 
to quantify, except for models Q and eventually by comparison with data. We 
shall assume that duality violations will be less than 10% for AF^^. 

To leading order in the heavy quark expansion, the long distance contri- 
butions to AFb^ are parameterized by the matrix elements of two dimension 
six operators 

Q = (iiSi)v-A(bjSj)v-A, (7) 
Qs = {biSi)s-pibjSj)s-p (8) 
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between a B„ and -B,, state. We write these matrix elements as 



(Q) 
{Qs) 



5 o o M'l 



(9) 
(10) 



where fs^ is the Bg decay constant. The 'bag' parameters B and Bs are 
defined such that B = Bs = 1 corresponds to factorization. B also appears in 
the mass difference, while Bs is specific to AFb^. 

The matrix elements of these operators are not independent of mb- Their 
mft-dependence could be extracted with the help of heavy quark effective the- 
ory. There seems to be no gain in doing so, since the number of independent 
nonperturbative parameters is not reduced even at leading order in 1 /mb and 
since we work to subleading order in l/m,b even more parameters would ap- 
pear. The matrix elements of the local AB = 2-operators should therefore be 
computed in 'full' QCD, for instance on the lattice. 

Including l/mfo-corrections, the width difference is found to be 



(72 2 



^z)K2 {Q) 



(11) 



+ (1 + 2z) {K, - K^) {Qs) + S^rn + Sr 

where z = ml /ml and 

= NcCl + 2(7i (72 K2 = Cl . 
The l/mfe-corrections are summarized in 

k/m = {1 + 2z)[Ki{-2{R,)-2{R2))+ K2{{Ro)-2{Ri)-2{R2)) 
12z^ 



(12) 



1-4:Z 



Ki {{R2) + 2(i?3)) + K2 {{R2) + 2(i?3)) 



(13) 



The operators Ri and Ri involve derivatives on quark fields or are proportional 
to the strange quark mass TOs, which we count as Aqcd- For instance. 



-Ri = —ibiSi)s-p{bjSj)s+p, 

Tflfy 

R2 = -\{bi^pD''si)v-A{bjSj)v-A- 



(14) 
(15) 
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Table 1: Dependence of a, b and c on the fe-quark mass (in GeV) and renormalization scale 
for fixed values of all other short-distance parameters. The last column gives {Ar/T)g^ for 
B = Bs = 1 (at given fj.), fs^ = 210 MeV. 



TOfc 


M 


a 


b 


c 


(Ar/r)B. 


4.8 


nib 


0.009 


0.211 


-0.065 


0.155 


4.6 


nib 


0.015 


0.239 


-0.096 


0.158 


5.0 


nib 


0.004 


0.187 


-0.039 


0.151 


4.8 


2mf, 


0.017 


0.181 


-0.058 


0.140 


4.8 


mb/2 


0.006 


0.251 


-0.076 


0.181 



The complete set can be found in Ref. Operators with gluon fields con- 
tribute only at order [KQCD/fnbY- Since the matrix elements of the i?^, Ri 
are l/mf,-suppressed compared to those of Q and Qs, we estimate them in the 
factorization approximation, assuming factorization at a scale of order mb (A 
smaller scale would be preferable, but would require us to calculate the anoma- 
lous dimension matrix.). Then all matrix elements can be expressed in terms 
of quark masses and the Bs mass and decay constant. No new nonperturbative 
parameters enter at order l/nib in this approximation. 

The term Srem denotes the contributions from Cabibbo-suppressed decay 
modes and pengiun operators. They can be estimatedEI to be below ±3% and 
about —5%, respectively, relative to the leading order contribution. We neglect 
this term in the following numerical analysis. 



3 Numerical estimate 



It is useful to separate the dependence on the long-distance parameters Jb^ 
B and Bs and write (Ar/r)s^ as 



AT 



aB + bBs + c 



ft 



210 MeV 



(16) 



where c incorporates the explicit 1/mb-corrections. In the numerical analysis, 
we express T as the theoretical value of the semileptonic width divided by the 
semileptonic branching ratio. The following parameters are kept fixed: mb — 
nic = 3.4 GeV, = 200 MeV, A^^^ = 200 MeV, Mb, = 5.37 GeV, B{Bs 
Xev) = 10.4%. Then a, b and c depend only on rrib and the renormalization 
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scale /i. For some values of nib and fi, the coefficients a, b, c are listed in Tab. ^ 
For a central choice of parameters, which we take as mt, = 4.8 GeV, fi — mf,, 
B = Bs = land /b, = 210 MeV, we obtain (Ar/r)s, = 0.220-0.065 = 0.155, 
where the leading term and the 1/TOfc-correction are separately quoted. We 
note that the V — A 'bag' parameter B has a very small coefficient and is 
practically negligible. The l/mf,-corrections are not small and decrease the 
prediction for AF^^ by about 30%. 

The largest theoretical uncertainties arise from the decay constant fs^ 
and the second 'bag' parameter Bs- In the large- A'c limit, one has Bs ~ 6/5, 
while estimating Bs by keeping the logarithmic dependence on nib (but not 
l/mfe-corrections as requiced here for consistency) and assuming factorization 
at the scale 1 GeV gives Q Bs — 0.88. Bs has never been studied by either 
QCD sum rules or lattice methods. In order to estimate the range of allowed 
AFb^ conservatively, we vary Bs = 1 ± 0.3, fs^ — (210 ± 30) MeV and obtain 

(^)^ =0.16t°-. (17) 

This estimate could be drastically improved with improved knowledge of Bs 
and fs, ■ 



4 Measuring AF^^ 

In principle, both F^ and Th can be measured by following the time-dependence 
of flavor-specific modesid, such as Bs — > Dglv, given by 

e-^«* -t-e"^^*. (18) 

In practice, this is a-tough measurement. Alternatively, since the average Bs 
lifetime is predicted □ to be equal to the Bd lifetime within 1%, it is sufficient 
to measure either or Th- 

The two-body decay Bs DfDj^ has a pure CP even final state and 
measures F^,. Since and do not decay into <j) as often as Ds, the 
(f>(f)X final state tags a i?s-enriched B meson sample, whose decay distribution 
informs us about F^. 

A cleaner channel is Bs J/4'4'i which has both CP even and CP odd 
contributions. These could be disentangled by studying the angular correla- 
tionsH. In-practice, this might not be necessary, as the CP even contribution is 
expected Q to be dominant by more than an order of magnitude. In any case, 
the inequality 

Fl > l/riBs -> J/V"^) (19) 
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holds. CDpi has fully reconstructed 58 J/ip(t> decays from run la+Ib 

and determined t{Bs J/ip(l>) = IMtoHl ± 0.05 ps. Together with r(Sd) = 
1.54 ± 0.04 ps, assuming equal average Bd and Bg lifetimes, this yields 



Ar 



> 0.3 ±0.4, (20) 

Bs 



which still fails to be significant. 

In the Tevatron run II, as well as at HERA-B, one expects 10'^ — 10^ 
reconstructed J/'00, which will give a precise measurement of AF ■ 



5 Implications of non-zero AF^^ 
5.1 CKM elements 

Once AFs^ is measured (possibly before AA/^^ is measured!), an alterna- 
tive route to obtain the mass difference could usc.|this measurement combined 
with the theoretical prediction for (AM/AF)^^ lillld. The decay constant Jb, 
drops out in this ratio, as well as the dependence on CKM elements, since 
\iycbVcs) / {VtsVtb)\'^ = 1 ± 0.03 by CKM unitarity However, the dependence 
on long-distance matrix elements does not cancel even at leading order in 1/ mf, 
and the prediction depends on the ratio of 'bag' parameters Bs/B, which is 
not very well-known presently. We obtain AF/AM = (5.6 ± 2.6) • 10^'^, where 
the largest error (±2.3) arises from varying Bs/B between 0.7 and 1.3. 

When lattice measurements yield an accurate value of Bs/B as well as 
control over the SU{3) flavor-symmetry breaking in Bf^, the above indirect 
determination of AMb^ iu conjunction with the measured mass difference in 
the Bd system provides an alternative way of determining the CKM ratio 
iVts/Vtdl, especially if the latter is around its largest currently allowed value. In 
contrast, the ratio T{B K*^)/T{B — > {g,uj}j) is best suited for extracting 
small |Vts/Vtd| ratios, provided the long distance effects can be sufficiently well 
understood. 



5.2 CP violation 

The existence of a non-zero AF^^ allows the observatiDiL,of mixing-induced 
CP asymmetries without tagging the initial Bg or Bs uxl. These measure- 
ments are difficult, but the gain in statistics, when tagging is obviated, makes 
them worthwhile to be considered. The mass difference drops out in the time 
dependence of untagged samples, which is given by 

A+(e-r^* ± e-^«*) + A_(e-^^* - e"^"*). (21) 
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A- carries CKM phase information even in the absence of direct CP violation. 

In combination with an analysis of angular distributions, a measurement 
of the CKM angle 7 from exclusive decays governed by the b ccs or 
6 — !■ cus transition can be consideredlij. 
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